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Percolating networks based on interparticle tunneling conduction are shown to yield a logarithmic 
divergent piezoresistive response close to the critical point as long as the electrical conductivity 
becomes nonuniversal. At the same time, the piezoresistivity or, equivalently, the conductivity 
anisotropy exponent A remains universal also when the conductive exponent is not, suggesting 
a purely geometric origin of A. We obtain these results by an exact solution of the piezoresistive 
problem on a Bethe lattice and by Monte Carlo calculations and finite-size scaling analysis on square 
lattices. We discuss our results in relation to the nature of transport for a variety of materials such 
as carbon-black-polymer composites and RuC>2-glass systems which show nonuniversal transport 
properties and coexistence between tunneling and percolating behaviors. 
PACS numbers: 64.60.Fr, 72.20.Fr, 72.60.+g 



I. INTRODUCTION 

Tunneling-percolating network models have been sug- 
gested to properly describe transport properties of disor- 
dered systems in which tunneling coexists with percola- 
tion behaviors. Examples of such systems are provided 
by various carbon-black-polymer composites,EJ for which, 
close to the metal-insulator transition, the conductivity 
a follows a percolation-like power law of the form: 



(1) 



where Co is a prefactor, p is the concentration of the 
conducting phase, p c is the percolation, critical concen- 
tration, and t is the critical exponent.cl A clear indica- 
tion that transport is dominated by tunneling between 
carbon- black grains is provided by a large strain and/or 
volume sensitivity of ctBu Similar situations are encouru 
tered also in organic conductors-polymer composites,u 
and thick-film resistors made of metal-oxide conduct™ 
grains inhomogeneously embedded in a glassy matrixE'l 
These latter compounds display quite large piezore- 
sistive effects and their use for sensor applications is 
widespread 

Modelling such systems in terms of tunneling- 
percolating networks resides on the assumption that the 
exponential decay exp(— 2r/£) of the tunneling amplitude 
permits to retain only nearest neighbouring conducting 
grains. Since the tunneling factor £ oc 1/W depends 
only on the potential barrier V separating the grains, a 
modification of the tunneling distance r induced by an 
applied strain e can lead to a large variation of the to- 
tal conductivity, and so to large values of the piezoresis- 
tive coefficient d\og(a~ 1 )/de. For example both carbon- 
black-polymer composites and Ru02-based thick-film re- 
sis||||s can display piezoresistive coefficients up to about 

Although percolating network theories predict that the 
conductivity exponent t should be universal and equal 



to t ~ 2.0 for three dimensional latticesja tunneling- 
percolating systems display important deviations from 
universality. For example, values of t up to about 
t ~ 6 have bffn reported for some kind of carbon-black 
compositespcl while t ~ 4Ji,has been measured for RuC>2- 
based thick-film resistors.113 

Nonuniversality has been theoretically proposed to 
arise from specific conducting-insulating phases distribu- 
tions such as in the Swiss-cheese model where spheri- 
cal insulating inclusions are introduced in a continuous 
conductors An alternative explanation has been put for- 
ward specifically for tunneling-percolating materials by 
arguing that, if the composite has a sufficiently wide dis- 
tribution of tunneling distances, the distribution of bond 
conductances between conductive grains has a power-law 
divergence£3 It is well known that speh anomalous dis- 
tribution can lead to nonuniversality.113 

In this paper we argue that the piezoresistive effect, i. 
e., the sensitivity of a upon applied strain, could be a 
decisive tool to investigate the origin of nonuniversality. 
We show that if nonuniversality of tunneling-percolating 
systems is given by a diverging distribution of tunneling 
conductances, then the piezoresistive response diverges 
logarithmically at the percolation thresold, while at the 
same time the piezoresistivity anisotropy remains univer- 
sal regardless of the universality breakdown of the con- 
ductivity exponent t. We obtain these results by an an- 
alytic solution of the piezoresistive problem on a Bethe 
lattice and by Monte Carlo calculations and finite-size 
scaling analysis on square lattices. 

This paper is organized as follows. In the next section 
we introduce the lattice model and the relevant quantities 
concerning the piezoresistive response. In section III we 
solve analytically the piezoresistive problem for a Bethe 
lattice while in Section IV we investigate numerically the 
square lattice. A discussion of our results in connection 
with real materials is presented in the last section where 
also the conclusions are drawn. 
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II. THE MODEL 

Let us consider a random resistor network where the 
bond conductivity distribution is given by: 

P(9)=ph(g) + (l-p)5(g), (2) 

where p is the fraction of bonds with finite conductiv- 
ity g with distribution h(g). In the following, we im- 
plicitly assume that the temperature is high enough to 
neglect other processes than intergrain tunneling (grain 
charging effects, Coulomb repulsion). This is certainly a 
good approximation for thick-film resistors which shows a 
very weak temperature dependence of a already at room 
temperature. Hence, within a tunncling-percolation 
framework, h(g) depends upon the distribution of tun- 
neling distances between nearest neighbouring grains. 
For narrow distributions, h(g) is peaked around g — 
go exp(— 2a/£) where a is the mean intergrain distance, £ 
is the tunneling factor, and go has units of a conductivity. 
We assume that the main dependence upon the mean in- 
tergain distance a is all contained in the argument of the 
exponential and for simplicity we set go = 1. In a first 
approximation, we can set 

h(g) = 6(g-e- 2a /Z), (3) 

as an extreme case for narrow distribution for tunneling 
distances. In this situation, close to the metal-insulator 
transition, transport follows Eq.(|l|) with universal expo- 
nent t depending only upon dimensionality and all de- 
pendence upon tunneling distance a is embodied in the 
prefactor ao which measures the average bond conduc- 
tivity. Hence, we expect that a change of the tunneling 
distance a induced by an applied external strain would 
affect only oq. 

Let us consider now the case in which TV spherical par- 
ticles of radius Ro are placed at random so that the tun- 
neling distance r between two neighbouring grains fluctu- 
ates. When the average distance a between two adjacent 
grains is such that a ^> R , Balberg argued that the 
salient feature of the resulting bond conductivity distri- 
bution is captured by a power-law divergence at small 
gE 

h(g) = (l-a)g- a , (4) 

where a = 1 — £/2a and the prefactor I — a assures the 
correct normalization .-of h(g) and p(g). As first shown 
by Kogut and StraleyEJ for values of a larger than some 
dimensionality dependent critical value a c , universality 
breaks down and the transport exponent t becomes (in- 
dependent: t — > t(a) > t. The important point here 
is that since a depends on the mean intergrain distance 
a, then the exponent t(a) can be affected by an applied 
external strain e and this would lead to an anomalous 
piezoresistive response. 

To investigate this issue, let us assume that the net- 
work is embedded in a homogeneous elastic medium 



and that the elastic coefficients of the network and the 
medium are equal. Moreover we set the bond directions 
parallel to the axes of a D-dimcnsional cubic lattice. In 
this situation, if we apply an uniaxial strain e along, for 
example, the rc-axis, then the tunneling distance for a 
bond directed along x changes to a — > a x = a(l + e), 
while the bonds directed along the other orthogonal axes 
remain unchanged: ai = a for i x. The strain- 
induced change of the tunneling distances leads therefore 
to anisotropic bond conductivity distributions p x (g) and 
Pi^x(g) — pig)- For an external electric field Ei applied 
along the z-axis and up to linear order in e, the result- 
ing network conductivities <7j are a x = <r(l — riie) and 
Ci^x — cr(l—T±e) where a is the unstrained conductivity 
and 

r„ = ^; ix = ^ (M*), (5) 

" de de 

are the longitudinal and transverse piezoresistive coef- 
ficients, respectively. These two quantities completely 
define the piezoresistive properties of the network. For 
example, the isotropic (or hydrostatic) piezoresistive re- 
sponse r obtained iay applying equal strain e along all 
bonds directions is:Eil 

r = ^- =r " + (Z/2 - 1)r± ' (6) 

where Z is the coordination number, while informations 
about the tortuosity the current has in flowing through 
the network are given by the piezoresistive anisotropy 
factor defined as: 

rii -r ± 

X= JL r ^, (7) 
1 II 

which measures the degree of macroscopic transport 
anisotropy.EJ Close to the-Dercolation thresold, \ displays 
a power-law of the form:La 

X-(P-Pc)\ (8) 

where the exponent A is the same of that govern- 
ing criticality of the conductivity anisotropy A = I — 
a y /(Tx of randomxeaistor networks with anisotropic bond 
conductances. Ej'E£nilJ In the present case, in fact, bond 
anisotropy is induced by the applied uniaxial strain and, 
since A is independent of the degree of bond anisotropy ,113 
A and x have the same critical behavior with the same 
exponent A. 

Having introduced the main quantities defining the 
piezoresistive properties of random resistor networks, let 
us now discuss qualitatively the effects of universality 
breakdown induced by the diverging bond conductances 
distribution function Eq.(|J). The anisotropy factor \ 
and its exponent A cannot be accounted in a simple way 
without explicitly solving the anisotropy problem, which 
is done in the following sections. However, the hydro- 
static piezoresistive response T can be simply evaluated 
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by noticing that, by definition, it is obtained by changing 
the tunneling distance a to a(l + e) for all bond direc- 
tions. Hence T can be readily found by differentiating 
Eq.(nl) with respect to e: 



T a < a c 

r + (1 -a) In (^)t'(a) a > a c 



(9) 



where To = — dlna^ 1 /de. For a < a c transport is uni- 
versal and the piezoresistive response is governed solely 
by the strain-dependence of <tq in Eq.(|l|) leading to the 
p-independent factor Tq. Of course, distribution func- 
tions for occupied bonds like Eq.(Q) would trivially lead 
to the same qualitative result, i.e., a piezoresistive re- 
sponse independent of the concentration of occupied 
bonds. Instead, as soon as transport becomes nonuniver- 
sal (a > a c ), the tunncling-percolation model of Balberg 
Eq.([|) predicts a logarithmic divergence of the hydro- 
static piezoresistive response. Note that such an anoma- 
lous behavior is expected also for distributions more com- 
plicated than Eq.(^) as long as their asymptotic behavior 
for g — > has a power-law divergence with exponent de- 
pending upon the tunneling distance. 

In the next sections we shall verify the correctness of 
Eq. (|) and calculate the p-dependence of \ by consider- 
ing two quite dinstinct cases: the Bethe lattice model, 
which is paradigmatic of high-dimensionality lattices foe 
which transport is governed by mean-field exponents,t£l 
and the two-dimensional square lattice. 



III. BETHE LATTICE 

The conductivity problem on a Bethe lattice, or Cay- 
ley tree, model has been considered and solved in Refs.[ 
poj M. As we show below, for this model the concen- 
tration dependence of the piezoresistive response in the 
critical regime can be obtained analytically for both dis- 
tributions Eqs.(|3[^). For simplicity in the following we 
shall consider a Bethe lattice with coordination number 
Z = 4 (p c = 1/{Z - 1) = 1/3)). According to Ref.[ 
p2| , the distribution Pi(a) of conductivities from an ar- 
bitrary bond directed along i = x,y to infinity and the 
current distributions Ji(a) induced by an applied electric 
field Ei satisfy the following coupled non-linear integral 
equations: 

P x (ai) = J dgp x (g) J da 2 da 3 da 4 P x {<J 2 )P y {(T 3 )P y {a i ) 
7(o-2 + 0-3 + 0-4) 
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g + (72 + 0-3 + (T 4 



(10) 



J x {<Ji)= E x aiP x (ai) + J dgp x (g) j da 2 da 3 da 4: 



gJx{cr 2 ) P y (<7 3 )Py (g4) 
g + 02 + 03 + (T4 
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5(02 +cr 3 + 04) 



g + (72 + (7 3 



(74 



The correspondin g eq uations for P y (a) and J v (cr) are ob- 
tained from Eqs . (|Ty ,|llj) by substituting x with y (and 
vicerversa), and the total conductivities Oi are obtained 
from: 



C-OO 

-2 / dzJ X ( Z )Py 

Jo 



d[P X ( Z )Py(z) 

dz 



(12) 



where Pi{z) and Ji(z) ap.the Laplace transforms of Pi(<r) 
and Ji(cr), respectively.E3 

Let us consider first the case of a narrow tunneling- 
distances distribution, Eqs.(||,||), for which, in the 
absence of applied strain, the conductivity is a — 
6exp(— 2a/£)(p — Pc) 1 for \p — _pJ-.<C 1 where t = 3 is 
the universal transport exponent Ed The effect of an ap- 
plied uniaxial strain, e 7^ 0, can be readily found by 
following Ref. [ ^2|. In the critical region, we find that 
the piezoresistive coefficients Eq.(||) reduce to: 



IK-L) 



1 + H^(P-Pc) 



(13) 



The above result captures the essential physics at 
the basis j-of the piezoresistive response of percolating 
networks. For p > p c , Tu > T± because the strain- 
sensitivity is stronger for sample conductivities measured 
along the direction of the applied strain. As p moves 
towards p c , the current currying paths become more 
and more tortuous decreasing therefore the macroscopic 
anisotropy induced by e until at p — p c the longitudinal 
and transverse piezoresistive coefficients become equal.Efl 
From Eq. (|l3|) , we find that the piezoresistive anisotropy 
X, Eq. ([?]), goes to zero as: 
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(p-Pc)- 



(14) 



Hence the anisotropy exponent for the Bethe lattice is 
A = 1 in accord with previous results.E3 

Now, we consider how the piezoresistive response of 
a Bethe lattice changes when the bond conductances 
have a power-law distribution as in Eqs. (|||). In do- 
ing so, we generalize the procedure described in Ref. |l3| 
to the anisotropic bond conductance case of Eqs . (|To)-|l 2[) 
and perform an expansion in powers of e. To obtain 
the piezoresistive coefficients, it is sufficient to keep only 
terms up to linear order in e. Here we stress that the re- 
sulting piezoresistive response for the Bethe lattice model 
depends crucially on the sign of a. For a < 0, corre- 
sponding to the quite unphysicaLrelation a < £/2, con- 
ductivity is universal with t = 3,Ej and the piezoresistive 
response can be easily shown to behave qualitatively as in 
the binary distribution case discussed above (finite values 
of r || at p = p c and A = 1). Instead when a > 0, the 
piezoresistive response changes qualitatively as discussed 
in the following. 

From Eq.(nG) we find that for a > the p — p c de- 
pendence of the Laplace transforms of P x (a) and P y (a) 
As: 
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P x{y) {z) =l-6f(u) + 61n5[g(w) + (-)5 gi (u;)}s, (15) 

where u) = 8 1 -" z and 5 = (p — p c )/Pc- The functions 
f{y)i 9{y) and gi{y) satisfy coupled integral equations, 
but their explicit expressions are not of interest here. By 
applying the same procedure to Eq.(|TT|) we find: 

■W*) = \s l+T ^E x {f'{^) - M</M + (-)6g 2 (co)]e}. 

(16) 

Finally, by substituting these results into Eq. (|l2|) we ob- 
tain 



ao — a\£ 



H5- 1 )[a 2 + (-)a 3 S]e}, (17) 



where ao, . . . , &3 are positive functions of jjl For e = 0, 
we obtain the result of Kogut and Straley:E3 <r x (y) = a ~ 
(p — PcY^ with t(a) = 3 + a/(l — a), indicating that 
transport is non-universal for a > a c = 0. In this regime, 
the piezoresistive coefficients are obtained from the terms 
of Eq.(|l7j) proportional to s: 



- II (-L) 



Oi +ln((5- 1 )[a 2 + (-)a 3 5] 
a 



In 



1 



(18) 



We have arrived therefore at the result that as long as 
the tunneling-distance distribution is such that transport 
becomes non-universal, then the piezoresistive response 
diverges logarithmically as p — ► p c . This must be con- 
trasted with the finite value of T|| and T± at p = p c , 
Eq . ( |l3|) , obtained from the simple binary distribution 
Eqs.(|2||3|). From Eq.(|l8|) we obtain also that the piezore- 
sistive anisotropy Xj Eq.(J^), reduces to 



X 



2a; 

02 



5 oc (p-pc)- 



Hence x goes to zero with universal exponent A = 
spectively of the logarithmic divergence of rum 
the universality breakdown of a. 



(19) 

1 irre- 
and of 



IV. SQUARE LATTICE 

The results of the last section have been derived for 
a Bethe lattice and are therefore in-general relevant for 
high-dimensionality (D > 6) latticesO To investigate the 
piezoresistive response for low dimensional networks, we 
consider now a square lattice. We employ the transfer- 
matrix method of Derrida and Vannimenus which per- 
mits to calculate exactly the conductivity of a long-strip 
of width N (along x) and length L~^> N (along y).c3 The 
distribution of the N x L conductors follows Eq. (§,§) and 
the longitudinal and transverse piezoresistive coefficients 
are obtained by setting a x = a(l + e), a y — a and a x = a, 
a y = a(l + e), respectively, with e = 0.001. For each 
value of a > at the percolation thresold p — p c = 1/2, 
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FIG. 1. Unstrained conductivity a (a) and piezoresistive 
anisotropy \ (b) calculated at p = p c — 1/2 as function of the 
width N of the strip for a = 0.05, 0.1, 0.15, 0.2, 0.3, ... 0.6. 
(c): solid circles: conductivity exponent i as a function of 
the distribution exponent a. Empty circles: results from Ref. 
p5| . Dashed curve is the theoretical estimate t(a) = 1.3158 
for a < a c = 0.24 and t(a) = 1/(1 — a) for a > a c . Solid 
squares: conductivity (or piezoresistivity) anisotropy expo- 
nent A. Solid line: average value A = 0.87. 



we perform the calculations at N = 8 (L = 1 x 10 8 ), 
N = 16 (L = 7.5 x 10 7 ), N = 32 (L = 2.5 x 10 7 ), 
N = 45 (L = 2 x 10 7 ), N = 64 (L = 1 x 10 7 ), N = 
82 (L = 1 x 10 7 ), and = 100 (L = 0.8 x 10 7 ). 

In Fig.^ and Fig.|^ we show our results for the un- 
strained conductivity a, the longitudinal and transverse 
piezoresistive coefficients r || , and the piezoresistive 
anisotropy X- From the a-vs- N data of FigJ^(a) we ex- 
tract the conductivity exponent t from the finite-size scal- 
ing relation a(N) = ciiV _t / ,y (l + GaViV) where v = | 
is the correlation length exponent O The resulting t- 
values (solid circles) are reported in Fig.|Tj(c) as func- 
tion of the power-law distribution exponent a of Eq. (|]) . 
For comparison, we report also the t values obtained in 
Ref. |25| by a different numerical method (empty circles). 
As expected, for a sufficiently large, the i-exponeiit in- 
creases well beyond its universal value t ~ 1.3,LJ'E3 sig- 
nalling breakdown of universality. According to Refs. 
|25| , pq , universality is lost for a > a c — 0.24 ± 0.08, 
and the conductivity exponent should follow the rela- 
tion t(a) = 1/(1 — a) reported in Fig.|l|(c) by the dashed 
curve. Our numerical results, as those of Ref. ^5], agree 
well with the theoretical expectations for large values of 
a but somehow overestimates the t- values for a ~ a c . 
This deviation from the theoretical expectations should 
be due to finite-size effects |aaialified by the presence of 
two competing fixed points.Prr 6 ! 

In contrast to the large a-dependence of <r(N), the 
data for different a values of the piezoresistive anisotropy 
X reported in Fig.|l|(b) collapse all in a single curve. 
Hence, the piezoresistivity (or conductivity) anisotropy 
exponent A [filled squares in Fig.^(c)] estimated from the 
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FIG. 2. (a) and (b): longitudinal and transverse piezore- 
sistive coefficients as a function of N for different values of 
the power-law distribution exponent a. (c): prefactor of 
Tikx) = a ||(±) ln(iV) as a function of a. 



finite-size scaling relation x(^0 — c^N^ x ^ u (l + c 4 /A) 
does not show appreciable variations over the entire range 
of a- values considered£3 This is in striking contrast with 
A = t — Pbi where /3 b — 0.48 is the fraction of conduct- 
ing bones in the backbone, conjectured to hold true for 
D = 2 in Ref. [H]. As for the Bethe lattice case, the 
universality of A also for a > a c suggests that this expo- 
nent depends only upon the geometry of the conducting 
cluster and is not influenced by t. 

The results for the piezoresistive coefficients T\\ and T± 
as functions of N are reported in Fig.|^(a) and Fig.|2|(b), 
respectively. With the exclusion of the large N values 
for a = 0.6, which we think are affected by too small 
values of L, Fii and T± follow approximatively a ln(iV) 
behavior. This signals a ln[(p — p c ) _1 ] divergence of the 
piezoresistive coefficients as p — > p c . How the logarithmic 
divergence depends upon a is studied in Fig.^](c) where 
we plot the prefactors au and a± of the finite-size scaling 



law: T 



l(J 



: a ||(_L) ln(iV)(l + •••). The correction-to- 
scaling terms • • • which best fitted the data where pro- 
portional to 1/N and to l/ln(iV). Clearly, aiim is a 
monotonous increasing function of a indicating that the 
logarithmic divergence of r || is stronger for a larger. 
In analogy with the results on the Bethe lattice, we would 
expect that oiim vanishes for a < a c , so that Tun) 
has a finite limit at p = p c . However the GS|i(j_) data of 
Fig.^(c) are small but nonzero even for a < a c = 0.24. 
We think that this is due to the same finite-size errors 
affecting the unstrained conductivity exponent, Fig.[I](c), 
which lead to a spurious a-dependence of a and conse- 
quently to a||(j_) 7^ 0. 



V. DISCUSSION AND CONCLUSIONS 

Let us discuss now the applicability of our theory to 
real materials. The thermal expansion effect on the re- 
sistance of nonuniversal [t ~ 3.0) carbon-black-polymer 




FIG. 3. Conductance G (a) and piezoresistance K (b) of 
Ru02-glass composites for various Ru02 volume concentra- 
tions x. Open symbols: Ref. Solid symbols: Ref. ^. 

composite reported in Ref. ^ provides a first indirect clue. 
From the resistance- us- volume data we have extracted a 
hydrostatic piezoresistive coefficient T of about 30 which 
keeps increasing as the volume is expanded by the tem- 
perature. If volume expansion effectively reduces the 
carbon-black concentration, we obtain a piezoresistive co- 
efficient enhancement as the concentration p moves down 
to its critical value. However, due to the uncontrolled ef- 
fect of polymer melting on the microstructure of the com- 
posite, we have been unable to single out any logarithmic 
divergence of the piezoresistivity response. 

A much clearer situation is found for Ru02-based 
thick-film resistors. In Fig||(a) we report conductance 
G measurements on several Ru02-glass composites fog 
different values of the metal volume concentration xBxi 
When ln(G/G ) 1/ * is plotted as a function of ln(a; — x c ), 
where Go is a prefactor and x c is the critical volume con- 
centration, the whole set of data collapses into a single 
straight line indicating a power law of the form: 



G~G {x-x c y 



(20) 



In Fig.||(a), the different symbols refer to different rela- 
tive particles sizes of the glass and RuC>2 and the values of 
x c and t which best fit the experimental data are reported 
in the inset. Clearly, all composites display values of t 
much higher than the three dimensional universal value 
i~2.oH 

The origin of such universality breakdown is investi- 
gated in Fig.||(b) where we plot the cc-dependence of the 
longitudinal piezoresistance coefficient K ^dlxiG^ 1 / de 
obtained by cantilever bar measurements. The mea- 
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sured K values fit reasonably well with a logarithmically 
divergent function of the form: 



K ~ K + B In 



1 



x c 



(21) 



represented in the (K — Kq) / B-vs-\n(x — x c ) plot of 
Fig.||(b) by the straight line. Since all the samples lie well 
within the critical region [Fig.||(a)], we expect that the 
piezoresistive anisotropy x is sufficiently small to regard 
if as a good approximation-pf the isotropic (hydrostatic) 
piezoresistance coefficient.!^ Hence, the data of Fig.|| are 
fully consistent with Eq.(^) for a > a c and provide a 
rather good example of the effect we have described in 
this paper. 

Unfortunately wc are not aware of reported measure- 
ments recording \ as a function of a;, so that the uni- 
versality of A claimed here even for a > a c cannot be 
verified. However the measurements reported in of Ref. 
p9| show that x decreases as the sheet resistance of com- 
mercial Ru02-based thick-film resistors increases. This 
is in qualitative accord with x ~ ( x ~ x c) X if higher re- 
sistance values are due to lower RuC»2 concentrations. 

In summary, we have shown by means of analytical 
and numerical results that when the tunneling exponent 
a = 1 — £/2a of the power-law distribution Eq.([|) is 
such that transport becomes non-universal, the piezore- 
sistive response changes drastically leading to a logarith- 
mic divergence of the piezoresistive coefficients as p — ► p c . 
In addition, we have demonstrated that despite of the 
universality breakdown of transport, the conductivity 
anisotropy exponent A remains universal. These features 
seem to be quite robust and calls for experimental verifi- 
cations on systems like carbon-black-polymer composites 
and thick-film resistors for which tunneling-percolation 
mechanism of transport have haaa .proposed and non- 
universality has been reported.EJoE3Ej. Earlier expert 
imental results on carbon-black-polymer. composites,cl 
and especially RuC>2 thick-film resistors J3'Q seem to indi- 
cate that indeed these systems are in the diverging tun- 
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